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Preface 


This  report  is  the  result  of  the  fall  of  the  Soviet  Empire  and  the 
realization  by  the  USAF  that  the  opportunity  to  obtain  excellent 
technology  and  also  support  Russian  scientists  was  in  US  interest.  In 
response  to  requests  to  initiate  promising  contacts,  Ted  Saivi  of  Phillips 
Laboratory,  acted  on  a  suggestion  from  Prof  Tony  Siegman  and  invited  Dr 
Yurii  Anan’ev  to  visit.  Prof  Siegman  had  met  a  number  of  prominent  Soviet 
scientist  in  St  Petersburg  in  1991  and  understood  Phillips  Laboratory 
interests.  These  three  scientists  are  united  by  their  interest  in  and 
contributions  to  laser  optics,  particularly  unstable  resonators.  Prof 
Siegman  wrote  an  introduction  to  Dr  Anan’ev  ‘s  book,  "Laser  Resonators 
and  the  Beam  Divergence  Problem,”  which  should  be  consulted  for 
historical  context. 

The  government  to  government  invitation  finally  succeeded,  due  in  part  to 
the  exceptional  enterprise  of  Ted  Saivi,  and  a  month  long  visit  to  the  US 
was  arranged.  Dr  Anan’ev  (he  writes  his  name  in  english  as  Ananiev)  and 
his  wife,  Tatania,  came  to  Albuquerque  and  Phillips  Laboratory  the  week 
of  January  11,  1993,  and  lectured  at  PL/LIDA,  the  IEEE  Lasers  and 
Electro-optics  Society  meeting,  and  at  the  University  of  New  Mexico.  Ted 
Saivi  was  host  for  this  unprecedented  visit. 

Following  his  return  to  St  Petersburg,  Russia,  Phillips  Laboratory  made 
funds  available  through  the  ILIR  (In-house  Laboratory  Independent 
Research)  program  for  a  small  research  effort.  An  appropriate  problem 
had  been  identified,  and  a  contract  with  Dr  Anan’ev  and  associates  at  the 
Laser  Technology  Center,  St  Petersburg  State  Technical  University,  was 
established.  This  report  documents  that  effort. 

The  extensive  contributions  of  Harro  Ackermann,  Greg  Dente,  and  Charles 
Moeller  to  this  project  are  gratefully  acknowledged.  Visit  support  was 
contributed  by  PL,  SPIE,  HDOS,  TRW,  Rofin  Sinar  (Seimans)  and  LEGS. 
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US  Authors’  note: 

Dr  Anan’ev’s  reports  were  sent  from  Russia  by  fax,  the  most  reasonable 
means  at  that  time.  Even  though  the  quality  of  the  facsimile  is  low,  to 
insure  accuracy  we  have  chosen  to  reproduce  the  reports  as  received,  with 
the  source  material  sent  to  Russia  as  an  Appendix,  and  a  short 
introduction  and  summary.  The  responsibility  for  this  additional  material 
rests  entirely  with  the  US  authors. 


1.  INTRODUCTION 


Unstable  Resonators  Driven  by  Spontaneous  Emission 

A  curious  phenomenon  was  discovered  in  unstable  resonators  in  diode 
lasers  in  1992  and  1993.  The  spectrally  dispersed  near  field  of  the 
output  of  an  unstable  resonator  diode  laser  pumped  just  below  lasing 
threshold  was  observed  using  a  grating  spectrometer.  This  is  possible 
because  diode  lasers  are  essentially  one  dimensional  emitters,  the 
transverse  direction  (growth  direction)  is  a  very  narrow  single  mode 
waveguide.  An  unstable  resonator  can  be  created  in  the  lateral  direction 
by  several  methods,  and  an  active  research  effort  has  been  underway. 

The  startling  data  that  were  taken  when  the  resonator  was  just  below 
threshold  displayed  a  herringbone-like  pattern  with  anomalous  curvature. 
The  model  at  each  frequency  seems  to  have  two  major  bumps,  which 
change  continuously  as  a  function  of  frequency.  The  modes  are  below 
threshold,  but  spontaneous  emission  at  each  frequency  is  preferentially 
amplified  in  the  modes  seen.  A  paper  was  published  in  Electronics 
Letters2  in  September  1993  that  showed  the  phenomenon.  Greg  Dente3  did 
several  types  of  analysis,  all  of  which  are  intriguing  but  are  difficult  to 
tie  together  in  a  coherent  picture.  He  gave  a  talk  at  the  January  1994 
Winter  Colloquium  on  Quantum  Electronics  in  Snowbird,  Utah  on  the 
subject.  Excerpts  from  the  paper  and  viewgraphs4  from  the  talk  are 
included  in  this  report. 


1  Modes  of  unstable  resonators  are  not  proper  mathematical  modes 
as  noted  by  Anan’ev  in  paragraph  2  of  his  first  (interim)  report.  These 
‘quasi-modes’  are  useful  concepts,  however,  and  will  be  referred  to  as 
modes  in  this  report. 

2  Bao  et  al..  Electronics  Letters,  23.,  pp1 597-1 598,  (Sep  93) 

3  Dente  and  Tilton,  PL-TR-94-1061,  pp  92-103. 

4  Unpublished 
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The  most  dramatic  attribute  of  the  data  was  the  direction  of  the 
curvature  of  the  patterns.  This  curvature  is  opposite  of  the  patterns 
observed  in  Fabry-Perot  lasers.  The  curvature  in  Fabry-Perot 
lasers  is  easily  understood  by  requiring  that  the  z  component  of  the 
propagation  constant,  kz,  remains  fixed  for  each  longitudinal  mode  (kzL/jc 
is  the  longitudinal  integer  mode  index),  and  arbitrary  ky  provides  the 
shape  of  the  curves.  In  the  farfield,  the  higher  spatial  frequencies  have 
higher  total  k  and  thus  shorter  wavelength.  The  curvature  reverses  for 
the  unstable  resonator  lasers. 

The  data  raise  fundamental  questions  about  the  existence  and  nature  of 
spontaneous-emission-driven  higher  loss  modes  in  a  running  resonator, 
and  about  the  mechanism  at  work  during  the  turn-on  of  an  unstable 
resonator.  Professor  Anan'ev  was  the  ideal  person  to  work  on  this 
problem.  He  is  a  pioneer  in  laser  unstable  resonators,  and  his  books  js  a 
classic  in  the  field.  He  had  also  been  somewhat  exiled  in  his  own  country 
for  the  last  two  decades,  having  been  'purged'  from  the  USSR  military 
research  programs.  Last  year,  one  visiting  Russian  graduate  student  said 
that  Anan'ev  was  a  'legendary  figure'  in  Russia  and  that  his  book  was  hard 
to  obtain  there.  We  felt  we  would  all  benefit  by  working  together  on  a 
problem  so  basic  to  laser  operation  which  might  also  shed  light  on  the 
fundamental  limitations  to  coherent  laser  performance. 

The  reports  provided  by  Professor  Anan'ev  are  difficult  to  read.  This  is 
partly  due  to  the  nature  of  the  subject  and  partly  due  to  language 
difficulties.  To  make  this  report  more  readable,  we  have  added  some 
background  material:  a  Preface,  this  first  chapter,  an  Appendix  containing 
Dente’s  Viewgraphs  that  were  sent  to  Russia,  and  a  summary  chapter. 

The  second  chapter  is  the  original  report  made  by  Professor  Anan'ev:  An 
interim  report  provided  in  October  1994.  In  March  1995,  another  report 
was  submitted,  with  a  computer  code  which  calculates  the  onset  of  lasing 
in  a  resonator,  (Chapter  3).  Final  comments.  Chapter  4,  were  provided  in 
July  1995.  The  final  chapter  is  a  summary  of  conclusions  on  the  topic. 


5  Yu  Anan’ev,  Laser  Resonators  and  the  Beam  Divergence  Problem, 
Adam  Hilger,  Bristol,  Philadelphia  and  New  York,  1992. 
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EXPLANATION  OF  BEHAVIOUR  OF  SEMICONDUCTOR  LASER  WITH 
UNSTABLE  RESONATOR  IN  SUB-THRESHOLD  REGIME 
(INTERIM  REPORT^ 


1,  Experimental  daia 


Tie  experimental  data  concerning  a  properties  of  diode  laser  behaviour  in 
sub-threshold  regime  and  having  been  informed  by  Phillips  Lab  are  as  follows. 

The  observations  have  uncovered  an  unusual  property  in  the  output  of 
unstable  resonators,  it  is  possible  to  spectrally  disperse  the  near  field  oui^ut  of 
an  unstable  resonator  diode  laser  by  using  a  grating  spectrometer.  Since 
diode  lasers  are  essentially  one  dimensional  emitt«*s,  the  transverse  direction 
(growth  direction)  isav^'naiTuWsiiigjle  mode  waveguide.  An  unstable 
resonator  can  be  created  in  tiie  lateral  direction  by  several  methods.  The 
interesting  laboratory  data  was  obtained  when  die  resonator  was  just  bdow 
direshold.  It  consisted  of  a  pattern  similar  to  a  herringbone.  The  mode  at  each 
fiequency  appears  to  have  two  major  bumps,  vriiich  change  continuously  as  a 
function  of  fi^cquenc)^.  The  modes  are  below  threshold,  but  spontaneous 
emission  at  each  frequency  is  preferentially  ainplified  in  the  observed  modes. 
Besides  this  information,  there  arc  many  photos  of  spectrally  resolved  near-field 
distribution  of  diode-laser  radiation  in  sub-threshold  regime. 


It  is  well-known  that  if  for  two  neighbouring  resonance  frequencies  of  a 
sj'stem,  CO,  ■»  ©5  -  and  ©2  »  iQ2  (the  timefactor  of  the  complex-valued 

amplitude  is  exp  (~i©t)),  the  following  inequality  is  valid 


2  1  ©I  -  ^2 1  /  (  of  +  ©5  )  <  1 


then  th^r  rCSviiSncc  overinp  (5^*  fll^ 

corresponding  evaluation  for  tiie  first  and  fiie  second  modes  of  a  plane  resonator 
vviu'i  uic  Muiie  Ovciiui  uiuicilSiOuS  ailu  SuHi  lu5565  of  iauiation  pcT  uOufils  tTonSlt 
as  is  the  case  of  diode-laser  unstable  resonator. 


Using  Ac  notation  of  {!],  section  3. 1,  we  have  2|a}-fflj|/(a,'+«»5)=' 

=  4| 5|  - 82I /  (A,  +  Aa  )  where  5^  are  phase  corrections  and.  A,^  are  the  relative 


losses  of  radiation  intensity  per  double  transit  For  Ae  plane  resonator  modes 
®  ^  (“  ■where  N  is  Fresnel  number  and  m  is  transverse  mdex. 


Ilierefoie  |  -  S|| » —  ;  for  conditions  resonator  is  filled  by  medium 

with  refractiYe  coefiBciciit  is,  N  r-  /  ?Xy  vAere  a  is  ths  half  v^idtih,  X  -  wave 
length  in  vacuum,  and  Lq  -  the  lengft  of  resonator.  In  our  case  X,=0.97  um,  a  - 
100  um,  Lo“500  um,  thus  N  «72  and  1 8J>  -  sj  |  «0.0I7.  Because  of  larger  width 


of  beams  in  unstable  resonators  the  phase  corrections  in  this  case  are  still 
smaller.  At  the  same  time  the  losses  of  radiation  per  double  transit  even  in 
condition  of  their  partially  compensation  by  gain  in  sub-threshold  regime  turn 
out  much  higher.  In  actual  fact,  threshold  gain  per  double  bounce  is  equal  here 
M/R'  =  7.8  where  R'  is  tfic  Fresnel  reflection  coefficient.  If  intensity  of  the 
excitation  is  lower  titan  threshold  one  only  at  5%,  the  gain  is  equal  7.8®”  *7.0, 
and  the  relative  losses  for  the  first  mode  of  unstable  resonator  are 
(7.8  -7.0)/  IS  *  1094;  for  the  second  one  fliey  are  dramatically  higher.  We  can 

/"IN  - —  — -.J- 

UAOA  wAivuuvii  )  io  wviwudijr  vcuiu>  <uiw  uiiouuiiv  iC*ovuicu%Ji  uiimo 


structure  cannot  be  observed. 

The  unstable  resonator  modes  cannot  be  used  also  for  series  expansion  of 
arbitrary  tight  Aetd;  dgenfunctions  of  unstable  resonator  do  not  make  up  a 
complete  orthononnalized  set,  so  tiiat  one  cannot  use  such  an  mansion  [2, 1]. 
Therefore  we  are  forced  to  give  up  attempts  of  use  an  infomation  about 
resonator  modes  and  transfer  to  methods  of  classical  optics. 


3.  Explanation  of  phenomena  experirngitallv  observed  finm  die  viewpoint  of 
classical  optics 


The  law  that  describes  tiie  spatially-spectral  pattern  is  similar  to  the  one  for 
fiinges  of  equal  thickness  [3].  To  evaluate  this  depeaidenc^  you  have  to  find  out 
tile  distance  L(x)  fiom  plane  mirror  to  tiie  convex  one.  This  distance  has  to  be 
calculated  along  tiie  wiy  parallel  to  tiie  axis  and  remote  fix>m  it  at  x  (see  Figure). 
It’s  easy  to  find  that 


L(x)  =  Lo  + 


2R 


(2) 


wh^e  radius  R  =  •*  2200  um  (M“  2.5  -  magnification  coefficient  of 

(M-1) 

resonator). 

The  interference  at  tiie  axis  takes  place  for  wavelength  X,  : 


U*n  =  m-J!- 


.(j  ‘II  —  III — ^ 

where  n  -  index  of  refioction,  q  -  the  order  of  interference. 

‘Pnv  'v  C\ 

A.W4.  MAW  WWAAAVW  •»  •  W 


(3) 
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L(x).n  =  mi!2>l^aiW 


and  dierefore 


'  R-m  2Rl_o 


•31 


(4) 


(5) 


One  cm  attach  tiie  index  *‘q”  to  this  parabola  Tbe  experimental  spatially- 

specttal  pattern  is  fee  sequence  of  predsely  feese  such  parabolas.  Their  spectral 
distance  at  the  axis  is 

(6) 

A  straight  Ime  corresponding  to  ai^'  ?,  intersects  the  sequence  of  parabolas  in  the 
points,  forming  fee  succession  of  intet^ence  peaks. 

Thus  the  shape  and  position  of  spots  csqjerimcntalty  observed  can  be  ogilained 
by  fee  idea  of  fringes  of  equal  thickness. 

4.  About  spectral  and  spatial  widfe  of  pattfm 

^  In  to  section  we  shaU  fescuss  fectors  defining  fee  widfe  and  number  of 
arcs  on  spectral-resolved  distributions  'nhidi  are  observed  ©cperimaitally.  We 
don’t  know  airy  woite  considering  interferometers  wife  medium  which  emits 
and  amplifies  a  li^t  in  sub-threshold  regime.  An  interference  pattern  must  be 
clear  if  an  elementary  wave  is  passing  at  least  two  times  ferou^  the  same  place 
and  has  approximately  fee  same  amplitude.  Let  us  use  this  assumption. 

uj  K^v^ruu  xpifcirui  widih  of  vxp^rirnffrtiul [/uU^rri, 

We  see  prunanly  feat  fee  spectra!  widfe  observed  is  several  times  smaller  than 
one  oflumincsccnce  line.  The  explanation  is  very  simple;  wc  saw  feat  fee  gain 
coefficient  drops  rapidly  as  frequency  moves  away  fiiom  fee  centre  of 
luminescOTce  line.  Concurrently  fee  wave  intensity  after  two  bounces  drops 
also,  and  int&fercnce  pattern  does  not  keep  existing. 
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b)  Spatial  width  of  interference  pattern. 

Let  us  consider  ffae  conditions  of  interference  pattern  existence  when  angle 
between  mirrors  is  a. 

We  can  consid^  diat  our  elementary  wave  moves  after  two  bounces  on  the 
same  tr^ectory  only  if  condition  w  X/ hn,  or 

h£X/an  (7) 

is  valid,  v^ere  oc^  is  diffiaction  an^e  and  h  is  die  size  along  die  wave  ftont.  On 
die  odier  hand,  in  order  to  avoid  significant  losses,  should  be  not  excessively 
big,  and  condition  h,  or 

b*i2LgX/n  (8) 

must  be  valid. 

We  can  easily  see  that  (7)  and  (8)  can  be  valid  simultaneously  if  a  <  ^yZL^n. 
For  magnitude  ag  of  the  half-width  of  zone  wdiere  angle  between  mirrors  of  our 
resonator  is  smaller  dian  a  we  find  bq  s  OgR  »  3Sum .  Taking  into  account  the 
Fou^  of  evaluation,  we  can  conclude  diat  this  value  corresponds  to  the 
experimental  data  quite  satisfactorily. 

5.  Conclusion 

llie  results  ot  me  above-proviOed  calculations  penmt  to  conclude  mat  me 
nhennmena  model  offered  bv  u.s  has  exnlaineld  oualitativelv  all  exnerimental 

r - - -  ^  B,  ^  B 

data. 
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ADDITIONAL  CONSIDERATIONS  ON  THE  BEHAVIOR  OF 
SEMICONDUCTOR  LASERS  WITH  UNSTABLE  RESONATORS 

I 

OPERATED  IN  THE  SUB-THRESHOLD  REGIME, 

AND  A  PRELIMINARY  DESCRIPTION  OF  A  COMPUTER  PROGRAM 
FOR  CALCULATING  THE  TPvANSITION  OF  AN  UNSTABLE-RESONATOR 
LASER  FROM  SUB-THRESHOLD  OPErIaTION  TO  lASlNG 

(DETAILED  PROGRAM  STATUS  REPOR  T) 


Jhis  leporl  presents  some  additional  considerations  concerning  the 
problem  of  radiation  coherence  in  suh-threshqld  regime,  which  was  not 

regime  where  ilie  lasing  starts  fronir  we  believed  it  appropriate  to  study  in 
more  detail  this  process,  and  to  develop  the  corresponding  computational 

I 

(wnh-nimji^  fn  with  unstable  resonators. 


1.  On  the  IVobiem  of  Radiation  Coh^ffi 


^ffsacc  in  Sub-Threshold  Regime 


'fhe  preceding  report  described  the  observed  pattern  of  radiation 
intensity  spectral  distribution  over  the  resonator  output  mirror  in  terms  of 
classical  optics;  in  more  than  one  respect  it;  is  similar  to  a  many-color 
picture  produced  by  flluminating  a  thin  film  |  with  white  light.  There  is, 
however,  a  substantial  difference  associated  wi^  the  significance  of  mutual 

coherence  of  the  primary  radiation  sources.  I 

I 

Indeed,  in  the  case  of  thin  films  we  deal  iwith  an  illumination  source 
located  on  one  side  of  two  fairly  weakly-reflecting  interfaces.  Therefore 
this  phenomenon  can  be  adequately  described  in  the  simplest 
approximation  of  two-beam  interference.  Note  that  for  an  interference 
pattern  to  exist,  the  temporal  coherence  length  should  be  at  least  of  the 
order  of  twice  the  distance  between  the  interfaces,  hi  our  case,  based  on 

I 

Ihe  known  gain-profile  width  in  a  semiconductor,  Av  -  400  cm"^,  we  come 
to  a  coherence  length  of  OTder  25  pm,  which  means  that  this  condition 
cannot  be  met  in  principle. 


iLine  3  reads  “dealt  with  in  preceding  interim  report.  Since  it  is 
the  subthreshold  ...” 
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i 

Neverthelesst  senuconductor  lasers  clearly  ezhibst  the  preseoce  of 
interfereBce.  It  is  a  consequence  of  two  factors: ' 

(i)  the  emitting  atoms  reside  between  the  resonator  miirors 
characterized  by  large  reflection  coedtidrats^  . 

(ii)  the  medium  has  the  prcqp>erty  of  the  gain. 

To  use  the  language  of  quantum  mechanics,  the  presence  of  the 
mirrors  increases  the  probability  of  emissiim  along  the  resonator  axis,  the 
wavelength  of  this  radiation  meeting  satisfactorily  the  resonance  condition. 
Viewed  from  the  standpoint  of  dassical  optics,  means  that  the  resonant 
wavelength  produces  constructive  interference  |  on  the  resonator  axis  and, 
hence,  the  highest  intensity  of  emitted  light  (for  equal  einiited  power  per 

atom).  This  is  the  reason  for  the  spatial  dependence  of  the  amplified* 

} 

radiation  spectrum  and  for  its  on-axia  monochTomatization.  This  pattern  of 
the  phenomenon  suggests  that  It  should  be  observed  only  in  the  near*axial 
region,  where  each  atom  of  the  medium  produces  a  series  of  imag^  in  the 
minors  making  up  a  linear  periodic  emitting  structure.  The  linearity  of  this 
structure  breaks  down  as  one  moves  away  from  the  axis,  and  this  is  what  is 
responsible  for  the  destruction  of  the  coDsinictive  interference.  It  should 
be  pointed  out  that  the  significance  of  the  equiyal^t  sources  in  this  series 
falls  off  with  increasing  nimibeT  of  the  mirror  reflections  involved  the 
faster,  the  smaller  is  the  reflection  coeffident.  It  may  be  suggested  that  as 
the  gain  increases  and  approaches  the  lasing  threshold,  the  effective 
length  of  this  emitter  series  should  increase,  so  that  interference  should  be 
observed  within  a  progressively  narrow  ^al  region?  also,  since  the 
spectral  selecUvily  ol  Uiis  eiiuller  set  increases  in  the  axial  dixectian,  the 
spectral  lines  should  become  narrower,  These  considerations  can 
apparently  be  used  as  a  basis  for  devirfqpment  of  a  more  campr^ensive 
ineory  oi  uie  uuserveu  paiiHrn. 

2.  Computer  Program  for  the  Calculaticm  of  Lasing  Onset  in  Lasers  with 
Unstable  Resonators 

Our  computer  program  draws  upon  the  unique  mathematical  model  of 
t3)e  process,  which  was  develc^ed  by  ns  earlier  ‘and  based  on  an  analysis  of 
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the  spatial  evolution  and  amplification  of  the  lirminescence  emitted  at  each 
puint  of  the  excited  active  medium  in  all  direc^ons.  The  i^atial  evolution 
is  calculated  witliin  the  geometrical-opticB  approximation  (for  unstable 
resonators  with  a  magnificaticm  coefficient  M  2  this  approximatioD  is 
known  to  guarantee  a  satisfactory  accuracy).  To  cut  down  computer  time, 
the  medium  parameters  and  the  values  of  total  radiation  density  over  the 
resonator  cross  sectional  area  were  assumed  to  be  the  same  as  on  the  axis; 
our  experience  in  computing  the  characteristics  lasers  with  telescopic 
resonators  shows  the  errors  due  to  this  assumption  to  be  insignificant. 
Although  this  model  was  constructed  for  lasers  with  shc^  inversion  times, 
for  which  it  is  particularly  useful,  its  universality  and  usefulness  cannot  be 
questioned.  I 

Let  us  present  the  main  relations  needed  here.  In  the  case  of  a 
cylindrical  active  medium  with  diameter  2^  length  L  equal  to  the 


resonator  miiTor  separation,  we  can  write  for  the  radiation  density , 

’  TJ, .  ( i  u-H'i 

pu-,/,'  =  x:  ^  ^  y  r';- " — r-  J 


X 


fit  r 


I-iere  z  is  the  cuordinaie  tuuug  the  jesonatar  aria  (with  the  origin  at  the 
concave  irarror),  f  is  the  iiiiie  leckuiied  huui  the  pimip  turn-on  (and,  herxce, 
from  the  appearance  of  spontaneous  radiation  do  lasing  firequency),  y  is  the 
volume  density  of  spontaneous  radiation  power,  By,  is  the  ray  matrix 
element  of  the  optical  system  reprQRpnting  a  section  of  the  resonator- 
oqui valent  optical  line  between  the  cross  sections  lying  at  the  coordinates 
z  ^nd  z  Tp  is  one  half  C3f  the  zone  diameter  anj  the  £  plane,  which  is  seen 
from  the  observation  point  at  z  and  was  calculated  by  means  of  this  matrix 
with  due  account  of  beam  obstructions,  and  is  the  gain  coefficient. 

This  equation,  for  whose  derivation  the  reader  can  bo  referred  to  [1], 
should  be  solved  together  with  the  equations  describing  the  change  in  the 
state  of  the  mediiwn  acted  upon  by  the  pumping  and  by  the  radiation  field 
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k  ; 

with  density  p  initiating  stimulated  emissicm.  The  state  erf  the  medium,  in 
its  luTiir  governs  lae  quantilies  y  and  kgn 

This  system  of  integrodiffer^tid  equations  lyi  the  power  W^of  the 
output  radiation  and  its  angular  distributiim; 

j 

H/(/)  =  \F[z\L,t)ek', 

r  i.  r 

0[fp,f)=  \  F[z',L,t)T\[z*,<^dzf\  \  F{2’,i.,f)  (iF  , 

•  I 

where  ©(v,  # )  is  the  fraction  of  radiation  powejr  ccHifined  to  a  fai-field 
solid  angle  (p , 

Tl(z',tp)  = 

Bfi  and  are  elements  of  above-mentioned  matrix  for  z  —  L. 

In  our  earlier  work,  we  develt^ed  the  eprrespooding  computational 
program.  It  was  written  in  the  presently-obsdlete  programming  language 
designed  for  the  BESM-6  computers  in  use  at  the  time.  While  in  order  to 
cut  down  computer  time,  it  was  subjected  ito  some  simplificatians  in 
addition  to  those  mentioned  above,  it  was  fonn^  to  be  fully  adequate  in  the 
interpretation  of  experimental  data  obtained  on  copper-vapor  lasers  jl]  and 
excimer  lasers  [2],  as  well  as  when  designing  lasers  of  these  types, 

pcifcnned  within  Grant  has  included  writing  computer 
programs  in  Fortran  and  Pascal,  which  are  much  better  in  all  respects.  Our 
calculations,  in  particular,  use  exact,  rather  than  length-averaged,  values  of 

y  and  k^,  the  resonator  length  is  divided  into  40  sections,,  and  not  10,  as 

i 

before,  and  so  on. 

The  program  consists  of  several  modules!  cenespending  to  individual 
elements  of  the  system  and  different  stages  ofjthe  solution.  This  permits  us 
to  modify  It  at  will  for  use  with  different  respnator  parameters,  pumping 
iecliuiqijes,  and  kinetic  modols  cf  the  n»divuT»- 


10 


5"  ■■  i  . 

Data  are  inpuUed  by  straigliuorvirard  fils  sdijtisg  Isy  mease  of  a  special 
Cd'tCyT  bailt  into  the  program  asrt  started  whenever  required.  Output  files 
receive  output  data;  a  subprogram  for  viewing;  graphs  on  the  monitor  is 
naturally  also  provided* 

Ihese  programs  are  presently  in  the  stage  of  check-out  using  the 

I 

familiar  copper-vapor  laser  model.  In  the  next  stagOf  the  check-out  will  be 
completedr  and  the  various  modificafions  of  the  programs  to  deal  with  the 

j 

geometry  characteristic  of  seimcanductor  lasers  will  be  developed,  after 
which  the  programs  win  be  submitfed  to  Philips  Lab. 

The  problems  associated  with  the  possible  role  of  a  number  of  effects 
whose  inclusion  would  make  the  programB  too  Ibulky  even  for  present-day 
computers  pose  severe  difficulties.  It  stiU  remaihs  unclear  how  smaff  is  the 
value  of  M  at  which  the  crossover  from  geometrical-opticB  to  diffraction 
approximation  becomes  necessary;  what  is  the!role  of  radiation  coherence 
effects  in  the  coiset  of  lasing;  to  what  extent  the  presence  of  a  low-O 
resonator  affects  the  spatial  distribution  of  nois4  radiation,  and  so  on.  Some 
or  inese  piublt^xiia  be  in  th;?  ctngc  cf  r\f 

calculations  with  experimental  data. 

i 

! 
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CONSIDERATIONS  ABOUT  SUMMARY 
“ASK  PATTERNS  &  MODE  FORMATION  IN 
UNSTABLE  RESONATORS”  (OF  G.DENTE  AT  AD^ 

AND  INFORMATION  ABOUT  PROCEEDING  PROGRESS  IN  CREATION  OF 
PROGRAM  OF  LASING  ONSET 
(INTERIM  REPORT) 

I  *  .CQPiniBnta  iinnn  rwnnrt.rASH  iiaHBTiia  *  morta  fAnnafin^ 

wwwiAKItt  TowmmtnTK"  tnt  G.Dnnf«»  jtf  »H 

A.,  PtHiMninaxy  rwinarts 

«)  Despite  of  general  appreciation  to  Siegman,  one  can  not  say  that  he 

known  much  earlier  than  publication  of  Siegman's  paper  (of  1965);  ho 
formulated  only  suppoaition  that  in  spite  of  conunon  fho 

unstable  resonator  can  turn  he  nsefnl. 

H)  RnsiiltR  of  calculation  in  aeometrical  approzinistios  coincido  with 
''®®^®“poudiag  formtilu  jjs  the  hoot  "Laser  Hesooators  and  the  Beaui 
Divergence  Problem”  of  Anan'ev  (p.  U4.  line  18  feom  top)  after 

replacing  L  — >  2Z  and  R  ~  —Rm- 

c)  General  view  of  experimental  ASE  patterns  (Spatio- Spectral  Patterns)  is 
in  accordance  with  considerations  provided  in  our  first  report.  Indeed, 
much  below  threshold  (current  of  180  ma  or  220  ma)  we  seo  phenomena 
of  classical  optics,  near  threshold  and  over  it  the  pictures  correspond  to 
principal  ro&onator  modo. 

^ TgflsneyativQ  ainplifior  propertios  of  unstable  r^Bonatori^ 

This  section  is  very  interesting,  hut  reading  is  aggravated  because  of  too 
short  tract.  For  example,  sense  of  "t"  is  not  explained;  one  does  not 
point  ont  whether  near  field  or  far  field  is  provided  on  all  figures,  and 
80  on. 

Method  of  Neumann  series  is  fully  appropriate  to  analyzing  the  problem.  We 
notice  only  that  more  obvious  deduction  of  formula 


2  line  A  a)  2  reads  “invented  the  unstable  resonator.  Really  unstable 
resonators  were  well  ...” 
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r  &  -1  iss  lac  KubKlilu&laa  laisiesis  id  F  IS  tSitf  liya: 


part  of  expreasion  F(x)  =  tS(x)’¥}!i*F  tlio  saino  right  part  of  the  same 
expression,  then  snbseqnont  analogous  suhstitntion,  and  so  on. 

U  is  not  dsar  why  the  r.unusiary  coxnprises  a  c^sideiaticn  of  case  of 
plane  resonstoT;  it  is  not  .-rlsar  too  why  on  the  last  Isss  of  p.l2  we  hsve  sign 
r&thoj  th&n  ***« 

The  numerical  rosuHs  ore  useful  and  clear,  but  scone  qualitative  laws 
was  published  before  in  small— known  (for  Western)  papers  of  kassiaii 
sdentiftts. 

Namslyt  it  was  shows  in  |!,  2J  that  converging  wave  of  iifistable 
resonator  possesses  maxiimnn  gain;  relevant  evaluation  was  made  in  the 
same  works.  Subsequent  discussion  on  this  theme  was  published  in  [3  — Sh 
The  answer  to  the  question  on  beginning  of  p.  13  is  provided  for 
arbitrary  resonators  in  work;«  of  UeTenlMirg  <{G— 8|  and  others).  He  solved  by 
standard  variation  method  namely  the  same  problcmi  which  is  formulated  on 
p.]4  of  Summary. 

Finally  we  notice  that  it  was  shown  in  some  paper  |9,  10]  that  all 
attempts  to  pass  converging  wave  through  unstable  xeBonator  lied  to 
ilraraalic  increase  In  value  of  output  angular  divergence.  It  can  be  explained 
in  such  manner.  If  we  want  to  achieve  a  high  energy-transformation 
etncieiicy,  our  amplitiGr  have  to  v/erk  in  regime  o!  deep  saturaiion.  In  this 


iroluiu  o»  gaiu  ou  sovurul  lust  puSsos  thtcugh  resensior  is  no  high,  and 


fluxes  of  energy  on  this  steps  of  converging  wave  Avnlntion  are  nearly  equal. 
At  the  sanrie  time  the  c<m verging  wave  has  the  diffraction  value  of  divergence 
only  on  sole  last  stop;  on  other  steps  it  has  smaller  area  of  cross  section  and 
therefore  the  bigger  divergence. 


2^NFQRMAT10N  AttOtJT  PROCnF.nWG  PROfiRKSS  IN  CUnATION  OP 
imonuAM  nn  T  AgiTjn  nivgFT 


We  have  created  a  program  for  calculation  of  laser  onset  in  the  case 
when  an  active  medium  takes  up  only  a  pari  of  resonator  length.  This  one 
will  be  provided  in  our  niaiil  ic-port.  At  present  we  are  isfeuming  about  some 
results  of  calculations  having  boon  done  by  this  program. 
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All  ilgiiroB  relate  to  the  foliowing  equations  and  values  of  lasor 
parameters; 

«  -  17  .-^  =  ^(N,-N,)+-a  . 


N,-N,  p,  .  N,  m,  0,,,  N, 

N, +Ng  + N,  =  N  =  const  •  • 

light  donsity  p(z,  t)  =  *If  (z',  z,  t)  dz'  , 

•>ei 

p  (*  >■•’)-  I  (a  ^|i«j  , 

ontpnt  power  W(t)  =  c»a*fl-  ~1  jfpfz',  z,  t)  dz'  , 
gain  coefficient  k  =  o(n,  -  N  j  , 

m 

power  of  spontanoouK  radiation  source  v=  hvNj/t ; 

at  last,  active  zone  radius  a =0.5  cm,  active  zone  length  40  cm, 

I'v  =  i5-  iO'*  cm  ’  (is  the  some  cases  2*  lO”  cm"*),  c=  10“**  cm*, 
active  medium  temperature  1600<N^  (its  value  effects  on  Initial  level 
population  with  boaring  in  mind  the  statistical  weight  1,  6  and  4  for  level 
1 .  2.  3  respectively),  t  =  770 •  10"®  s  ,  T*  =  10"*  s.  X *=  0.51  p,  resonator 
m&gniScatim!  coafficisnt  M=30,  resonator  length  ft’om  40  to  00  sm. 

On  fig.  1  are  provided  curves  for  cases  of  two  difloront  programs  of 
calculation  -  with  averaging  gain  coefficient  over  the  active  medhiro  len0h 
and  without  it  .  One  can  easily  see  that  difference  exJsta  bet  is  not  too  large. 

Fig.  2  and  fig.  3  correspond  to  csror  of  N  =-  25*10'*  cm"’  and 
N  =  2’  10''’  cm  ’  rospoctiifely.  One  can  see  the  incToase  in  atom  density  leads 
to  dramatic  increase  in  depth  of  time  modulation.  Fig.  4  and  5  manifest  the 
existence  of  the  same  law  in  the  case  when  active  medium  ffils  only  the  part 
of  resonator  length. 

Fig.  6  belongs  to  case  of  intermediate  resonator  length;  finally,  on  fig.  7 
IS  represented  temporal  spectrum  of  Inser  pulse  in  one  out  of  rouimes. 
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Figuro  1.  Laser  pulse  calculated  with 
averaging  over  the  active  medium  length  (2) 
and  without  it  (1). 
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Figure  2.  Laser  power  as  a  function  of 


Figure  3.  Laser  power  as  a  function  of 


time  for  *  40  cm. 


time  for  Ljc5«=40  cm,  N'^2xl0i5cm-3 
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Figure  6.  Laser  power  as  a  function  of  Figure  7.  Temporal  spectrum  of  laser 
time  for  Lres=50  cm  P''^s®=  ^  "  Lrcs’=40cm,  2  -  Lret-SOcm, 

3  —  Lf(,6*=*60cjii. 


program  of  calculation  for  kinetics  of  lasers  with 

TELESCOPIC  RESONATORS 

1.  INP.EXE  IS  A  FILE  FOR  INITIATION  OF  THE  PROGRAM 

PROGRAM  INP; 

{$M  4000,0000,10000} 

USES  CRT  ,  DOS  ; 

CONST 
RP  =  ’  - 
VAR 

R  TEXT; 

VARFLG:  BOOLEAN- 

BEGIN 

CLRSCR; 

wSSlilr  CALCULATION  FOR  KINETICS  OF  lASERS’l; 

wSSfi^r  p^'^^SCOPIC  RESONATORS-); 
wtaTELr>i(  Press  any  keW 

READKEY;  ' 

{$!-) 

ASSIGN(F,  TNPt.DAT); 

RESET(F); 

IH  lORESULTO  OTHEN  BEGIN 
^|^^(RP,  'HLE  ’INPLDA-r  NOT  FOUND.’); 

HALT 

END; 

{$1-1-} 

READLNIF.  M); 

R&ADLN(F.  L); 

RE.ADLN(F,  .A); 

RE.ADLN(F,  WL); 

CLOSE(F); 

CLRSCR; 

WRITELNIRP,  ’  PARAMETERS  OF  CAVITY  :V 

•  ^SONATOR  MAGNIHCATION  M  =',  M  :3:1); 

RESONATOR  length  LfCMl  =’  L -S-n- 
WRrrELN(RP, '  TUBE  R.ADtUS  JcMl  = ' '  A  -SM 1- 

WrSfHj/pd  •  WAVELENGTH  (MICRO:  'i  =’,  WL  :1:3)- 

SIS?  '  WANT  TO  CHANGE  DATA  ,  V  •.ES,N-NO)r-j; 

WHILE  NOT  KEYPRESSED  DO- 

CASE  RE.ADKEY  OF 

■Y*,  *?-,  CHR(13):  BEGIN  ‘ 

SWAPVECTORS; 

EXECCED-EXE*.  'INPl.D.Ar); 

SWAPVECrORS; 

IF  DOSERROR  <>  0  THEN  BEGIN 
WRITELN(RP,  ’UNABLE  TO  EDIT  FILE’1- 
RE.A>DLN;  '' 

ILALT(1) 

END 

END;  W : 

END; 

CLRSCR; 

{$1} 

ASSIGN(F,  ’INPa-DAT): 

RESET(F); 

IF^ORESULT  <>  0  THEN  BEGIN 
VRrTEU-4(RP,  'FILE  '’INPS.DAT’  NOT  FOUND.'); 
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READLN; 

H.M.T 

hNi-; 

iii-rj 

READLN(F,  SIGMA); 

RII^LN(F,  N): 

READLN(F,  TAU); 

CLOSE(F); 

CLRSCR; 

WTONTO  •  PARAMETERS  OF  MEDIUM;’): 

“O®  SECTION  (lE-U’CM2).’,SIGM.A;4a|: 
„£U/CM3)  N  ;42|: 

’  UPPER  LEV^  LIFETIME  (NS)  TAU  ;4:2)- 

WRimilf^' '  CHANGE  D^A  (Y-YES,N-NO)  ?  *); 

WHILE  NOT  KEYPRESSED  DO 
C^SE  READKEY  OF 
"T,  y,  CHR(t3):  BEGIN' 

SWAPVECTORS; 

EXECCED-EXE*,  'INP2.DAT3; 

SWAPVECTORS: 

IF  DOSERROR  <>  OTHEN  BEGIN 
WRnELN{RP,  'UNABLE  TO  EDIT  FILE*)- 
READLN; 

HALT(l) 

END 

END; 

END; 

PUMHNG  ( 1/TH)  VARIANTS  ARE  SUGGESTED 

rr  andth  in  ns)-): 

WRITELNIRP^'INPUT  NUMBER  OF  VARIANT)- 
READLN(NV),-  '' 

IF  NV=  1  THEN  BEGIN 
W'RITELN(RP,  ’  Tt  =  ?  NS')- 
RE.^D(Tl); 

WRITELN(RP.  '  T2=  ?  NS'); 

READ(T2): 

.‘\SSIGN(F,  'NAK.DAT); 

RESET(F); 

REWriteif): 

WRITELN(F,  NV); 

WRrrELN(F,  Tl:4:4); 

WRITELNIF,  T2:4:4); 

CLOSE{F); 

END; 

IF  NV=2  "rHEN  BEGIN 
WRITELN(RP,  '  T1  =  ?'); 

RHADfri):  ih 

ASSIGN  (F,  ’NAK.DAT); 

RESETfF); 

R'EV/RrrE(F); 

WRITELN(F,  NV);  • 

WRITELN(F,  Tl:4:4): 

CLOSE(F); 

END; 

DT:=L/15; 

WRITELNC  TIME  OF  DOUBLE  CAVITY  P.ASS  =’  DT:2-3  ’  NS')- 
WRITELNC  input  OPER.ATING  P.ARAMETERS  OF  THE  PROGR-AM'!- 
WRITELNI’  number  OF  STEPS  ON  THE  CAVITY  LENGTH  (NL<=40)‘); 
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OJ  to 


READLN(N2]; 

NP:=  ROUND(400/2/NZ); 

'  nui^ber  of  double  cwitv-  d^s  -v 
\VRrrELN(’  NP<  =  ’,NP;3)-  ‘  ^ 

READ(NP); 

ASSIONfF,  ’N.DAT)- 
RESETfF); 

REWRrrE(F); 

WRrTELN(F,  NZ): 

WIUTELNfF,  NP); 

CLOSE(F); 

SWAPVECTORS; 

KJCEC(”rELESC.EXE',“): 

SWAPVECTORS; 

END. 

2.  input  parameters  are  located  in  HLES  -.DAT:  rNPl.INR2.N,NAK 
file  INP1.DAT: 

30.0  -  RESONATOR  M,AGNrFICATION 

40.0  -  RESONATOR  LENGTH  (CMl 

1.0  -  TUBE  RADIUS  (CM) 

0.51  -  RADLATION  WAVELENGTH  (MICRONS) 

FILE  INP2J3AT: 

1.0  -  STIMULATED  EMISSION  CROSS  SECTION  (1  E-U*rV5i 

20  -  MEDIUM  DENSITY  (LEU/CM3J 

770.  -  UPPER  LEVEL  UFETTME  (NS) 

HLE  N.DAT 

40  .  -  number  OF  STEPS  ON  CAVITY  LENGTH 
5  -  NUMBER  OF  DOUBLE  CAVITY  P.ASS 

HLE  NAKJJAT 

1  -  NUMBER  OF  PUMPING  VARLANT 

10.  -  PUMPING  RATE  IN  NS 

10.  -  -  PUMPING  INCREASE  RATE  IN  NS 

3.  RESULTS 

FILE  W{T).DAT 

3  SSS ;  X 

HLE  PHI(T1.DAT 

O^UT  POWER  FRACTION  WITHIN  THE  DiFRACnON  .AI4GLE  AS  A  FUNCHON  OF 

1  COLUMN  -  TIME  IN  NS 

2  ^OLUMN  -  TIME  IN  UNHS  OF  DOUBLE  CAVITY  OASS  TTMF 

3  COLUMN  .  OmPOT  POWER  FRACTION  WirS  TO^cnON  ANGIE 

file  w_phi.dat 

WTITHN  TOE  D^Sh^^GLE:*^^^  ENERGY  (mj)  IN  A  PULSE  AND  ITS  FRACTTC 
I  COLUMN  -  TIME  IN  NS 
COLUMN  -  PULSE  ENERGY  IN  tnj 

COLUMN  -  PULSE  ENERGi'  FR.ACTION  WITHIN  THE  DIFRACTION  ANGLE 
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KtALhUNCllON  VN(N Va'l/r2,Tj 
C  SUBROUTINE  OF  PUMPING  RATE  CALCULATION 
VN=0. 

IF(NV.EO:2)  VN=  L/Tl 
IF(NV.EQ.l)  VN  =  t./Tr (l.-EXP(-T/T2)) 

RETURN 

END 

Rid^  FUNCTION  RO(rZ) 

C  SUBROUTINE  OF  CALCULATION  FOR  AN  EFFECTIVE  MEDIUM  RADIUS  IN  A 
C  POINT  Z 

RE.^LM^ 

COMMON  /RESON/  AL.M,NZ,DZ;DT 
Z0=  (ABS(IZ)/Z/NZ)-INT(ABS(IZ)/2yNZ) 

RO=A 

IFfZ0.Gr^5AND.Z0.LT.l.)  RO=:A/M*(Z-M +2*(M-1.)*Z0) 

RETURN 

Em 

REAL  FUNCTION  DETB(IZ0,I4N) 

C  SUBROUTINE  FOR  CALCULATION  OF  A  BEAM  MATRIX  ELEMENT  B  FOR 
C  TRANSITION  FROM  Z  TO  ZO 
REAL  ML 

COMMON  /RESON/  AL.M,NZ,DZ;DT 

IF*(IZO.EO.IZ)  GOTO  1 

U=1 

IF(IZIXIZO)  U=2 
GOTO  (ll^l)rU 

I  DETB=0. 

N=0 

RETURN 

C  0=<Z0=:<L  Z=Z2>Z0 

II  F^iz/2Jrrz 
N=INT(F-.00001) 

T=0. 

IF{(N*ALOG(M)).LT.(10/ALOG(10.)))  T=ty(M**{2*N)) 

IF(N.EO.O)  T=L 

ir(F.GT.NAND.F.LE.(N  +  .5))  DETB=DZ*{NZ-IZO*f  (IZ-2/N*NZ)T+  NZ/(M 
P  -LJ^tl.-MT)) 

•  IF(F.Gr.(N-f  ,5)AND.F.LE.(N+  1.))  DErB=DZ‘(IZ-2.*N‘NZ-IZ0+  (M-L)/N 
P  Z"(NZ-IZ0)*(rZ-(2,‘N  +  10*NZ)  +  NZ/(M-l.)*(L-T)) 

RETURN 

C  0=<Z(K=L  Z=Z1<Z0 
21  F=-I2yZ/NZ 

N  =  INT(F+ .49999) 

T=0. 

IF((N^ALOG(M)).LT.(lO,‘ALOG(10.)))  T=  l./(M"n2"N)) 

IF(N.EQ.O)T=l. 

ir(r.GT.{N-.5) AND.F.LE.N)  DETB=DZ*(IZOT-(IZ+  (2.*N.1.)*NZ)  +  NZ/( 

P  M-l.)*(l.-'rMl) 

IF(F.GT.NAND.F.LE.(N  +  .5))  DETB=DZ*(NZ/(M-1.)*M*(1.-T)+T*(IZ0-(I 
P  Z+Z*N*NZ)  +  (M-1.)/M*IZ0*(IZ+Z*N*NZ)/NZ)) 

RETURN 

END 

REAL  FUNCTION  DETD(IZOJZN) 

C  SUBROUTINE  FOR  CALCUL.ATION  OF  A  MATRIX  ELEMENT  D  FOR 

C  TRANSITION  FROM  ZTO  ZO 
REALML 

COMMON  /RESON/  AL,M,NZDZDT 
IF-(IZO.EO.IZ)  GOTO  1 
J-  1 

IF(IZLT.IZO]  J=2 
GOTO  (1 1.21) J 
t  DETD=1. 

N  =  0 


22 


l  iNUc 
END  DO 

lF(W.NE.O.)  WPI=WPr/W 
IF(W.EQ.O.)  WPI=0. 

DO  jz==rr4 
IZl-NZ+JZ 

ROl = ROl  +  DCTF(NZ,IZl,rr)*DZ 
END  DO 
RO=ROl 
IP-(WXT.O.)  W=sO. 

IKfWPLLT.O.)  WPI=0. 
wRiTE(i,*)rrDT,rr/2./N2;w 
WRITE(2,")rrDTJT/2yNZ,WPI 
swsssw+wrD’ri^-a 

SP =SP+ WPPW^DT  l,E-3 
iF(swj'iE.o.)  WRrrE(3,*^)rrDT^,sp/sw 
RETURN 
END 

REAL  FUNCTION  DETFCIZ^ZlrTT) 

C  SUBROUnNE  FOR  CALCULATION  OF  RADIATION  FLOW  FROM  A  LAYER  DZ 
C  IN  APOINTZITOAPOLNTZ 
REAL  MU- 

COMMON  /RESON/  ALfM,N2;D2;DT 
COMMON  /MEDIA/  PHIO,SIGMAGaCOTAaSW,SP 
COMMON  /EVl/  Yl(Ch40.0:400) 

COMMON  /EV2/  72(0:40,0:400) 

BB=ALOG(ABS(DETB(IZ,IZLN))) 

DD=N"ALOG(M) 

G=:GAIN(IZJZLrr)  ‘ 

RR=ALOG(RO(IZl)) 

i=^rr-ABS(iz-izi) 

IF(LLT.0)  1=0 

F=  (ABS(IZl)/2./NZ)-INT(ABS(IZl)/2./N23 
lH(F.GE.0AND.FXE-5)  J=INT{P2*NZ) 

IH{F-GE-.5ANDPI-E.L)  J=INT(2*NZr(l.-F)) 

GAMMA=Y2(J.I) 

.  GG=ALOGCTAU) 

5EIF=GAMMA*EXP(G-GG-2-*(DD-RR+  BB))/4./30. 

RETURN 

END 

REAL  FUNCTION  GAIN(IZ,IZLrr) 

C  SUBROUTINE  FOR  CALCULATION  OF  RADIATION  O'JN  AT  PASS  FROM  A  POINT 
C  Z1  TO  A  POINT  Z 
REAL  MX 

COMMON  /RESON/  AL,M.NZ,DZ,DT 
COMMON  /MEDIA/  PHI0.SIGMAG0,C0,TAU,SW,SP 
COMMON  /EVl/  71(0:40,0:400) 

COMMON  /EV2/  72(0:40.0:400)  IP=rT-ABS(IZ-IZl) 

IF(IPXT.0)  IP=0 
S=0. 

10=1 

IF(IZl.GT.IZ)  IO=-l 

IZ2=IZl 

Doi=rp.rr-i 

F= (ABS(IZ2)/2./NZ).[NT(ABS(IZ2)/2./rte) 

IF(F.GE^0AND.F.LE-5)  J=INT(F*2*NZ) 

IF(F.GE..5AND.F.LE.L)  J=INT{2‘Nr  (l.-F))  S  =  S  + Y2(J.I)-Yl(J.I) 

IZ2^IZ2H“IO 
END  bo 

GAIN  =  S‘DZ"SIGM  A 
IF{nLEO.IZl)  GAIN  =  0. 

RETURN. 

END 
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NN=INT(200/NZ) 
rF(NP.GT.NN)  NP=NN 
WltlTEC/)’  STEP  IS  CHANGED' 

WRITE! V)'  N=',NZ 
II'(NZX£.40)  GOTO  999 
998  CLOSE(l) 

CLOSE(2) 

CLOSE(3) 

STOP 

77  IF(INT(IT/NZ/2)*2‘NZ£Q.IT)  WRITE(*,“)'  NP='JNT(rT/2/NZ) 

END  DO 
Cl-OSEd) 

CLOSE(2)  '' 

CLOSE(3) 

END 

REAL  FUNCTION  RO(IZTr) 

C  SUBROUTINE  OF  CALCUIATION  FOR  OUTPUT  POWER  AND  SATURATING 
C  INTENSITY  IN  A  POINT  (Z,T) 

REALN^i 

COMMON  /RESON/  Ai,M,N2.DZ,DT 
COMMON  /MEDIA/  PHIO,SIGMA,GO,CO,TAU,SWSP 
COMMON  /EVl/  Y1  (0:40,0:400) 

COMMON  /EV2/  Y2(0:40,0:400) 

ROl  =  Y2(IZTr)/30./TAU 
IFCIZXO.O)  GOTO  I 
IFCELEOJ'IZ)  GOTO  2 
DO  JZ=-IT,IT 

m=iz+jz 

IF(JZ.NE.0)  ROl  =R01  +  DETF(IZ,IZl,nTDZ 

END  DO 

RO=R01 

RETURN 

1  DO  JZ=-ITjr 
IZt=JZ 

IF(JZ.NE.0)  ROl=ROt  +  DETF(0,IZl,rD*DZ 
F^in  DO 
•  R0=R01 

RETURN 

2  W=ROl*3.14159*A*A*(l.-l./M/M)/2.*G0 
WPI=0. 

DOJZ=-IT,-l 
IZl=NZ+JZ 
CC= DETF(NZ,IZl,IT)*DZ 
ROt  =  R01  +  CC 

W=W+ CC*(l.-iyM/M)‘3.14 159*A*A*G0 
D=ABS(DETD(NZ,IZl,N)) 

B=ABS(DETB(NZJZ1,N)) 

IF  (BXO.O.)  GOTO  12 
IF  (D,EO.O.)  GOTO  11 
QF=0. 

IF((N‘ALOG(M)).LT.(8.*ALOG(10.)))  OF=iy(M**N) 
O2=2.*R0(IZl)*OF+2.*A*ABS(D) 

Q1  =  PHI0*B  ;,)i. 

0=1. 

IF(ABS(Q1).LTABS{02))  0=ABS((01/02))**2 
WPI = WPI  +  CC'O*(l.-l./M/M)*3.l4159'A*A‘G0 
GOTO  12 

11  Ol  =  N*ALOG(M)  +  ALOG(PHI0‘B) 

02 = ALOG(2.*RO(IZ1) ) 

0=0. 

IF(01.GE.02l  0=  1. 
rF!f.j.NE.l.)  6=EXP(01-Q2) 

WPI  =  WPI  +  CC'Q'd.-  l./M/M)'0'3.141 59*A*.VG0 
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VARIANT  OF  THE  PROGRAM  TELESC.FOR  TAKES  INTO  ACCUUN  i  uis,  i  mou  iiwin  v^r 
INVERSION  OVER  THE  CAVITY  LENGTH 

PROGR.^M  TELESC 
ClriARACTER*44  ANSWER 
RC.AL  M,L 

COMMON  /RESON/  ACM,NZJ>Z,DT 
COMMON  /MEDIA/  PHIO^IGMAGO.COTAU^,SP 
COMMON  /EVl/Yl  (0:40,0:400) 

COMMON  /EV2/  Y2(0:40,0:400) 

1  F0RMAT{F8.1A44) 

OPnN(UNrr=  l,FILE=lNPlJ>At') 

READ(14)MANSWER 
READ(l,l)UANSWER 
RaAD(14)AANSWER 
READ(  14  )WLANSWER 
CLOSE(l) 

OPEN(UNrr=  13LE=’INP2J>Ar) 

READ(14)SIGMAANSWER 
READ(14)C0hANSWER 
READ(  14)TAUANSWER 
CLOSE(l) 

OPEN(UNIT=  l,FILE=*NA£CDAT) 

RE.M:>(1/)NV 

READ(1,*)T1 

'12=0. 

IHNV£0.1)  READ(14)T2 
CL.OSE(l) 

OPEN(UNrr=  LFTTJE = 'N.DAT) 

READ(1/)NZ 

READ(V)NP 

CLOSE(l) 

S;CMA=SlGMA*C0 
PHr0=12E-4*WL/A 
G0=6.626l76*90./WL*CO 
999  DO  IX=0,r4Z 
-  Y1(IX,01=0. 

Y2(IX,0)=0. 

END  DO 

DZ=L/NZ 

DT=DZ/30. 

oPEN{UNrr=  uhle='w(T).dat) 

OPEN(UNrr=2,nLE = 'phi  (T).DAr ) 

OPEN{UNrr=  3,FILE = 'W_PHII>AT) 

sw=o. 

SP=0. 

DO  ET=1,NT 

T0=(IT-l)*DT 

rs=iT«i 

rND=0 

DOIX=0,NZ 

XI  =  l,-Yl(rX,IS)-Y2(IX,IS)  .. 

IF(X14T.0.)  Xl=0. 

IFCYl(rX(IS)XT.0..OR.Yl(IX,rS).GT4.)  Yl(IXIS)  =0. 

IF(Y2(IX.IS)  JLT.O..OR.Y2(rX4S).GT.  1.)  Y2{IX,IS) = 0. 

DEIU  =30/SIGrvLArROaX,IS)*(Y2(i:aS)-Yl(rX,IS)) + Y2(IX,IS)/TAU 

DER2=(Xl-Y2(IX,rS))*VN{NVTl,T2T0)-DERl 

Yl(IX,m  =Yl  (IX,IS)  +  DERrDT 

Y2(rx,rn  =  Y2(IX,IS)  +  DER2-DT 

IF((ABSOri(IX,IT))  +  ABS(Y2(IX.rr))).GT.l.)  IMD  =  I 

END  DO 

IF(IND.EQ.0)  GOTO  77  rF{SW.GT.l.)  GOTO  998 

nz=nz+i6 
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RETURN 

C  0<=Z0=<L  2=Z2>Z0 
It  F=I2y2./NZ 

N=INT(F-.00001) 

rF(F.GT.NJ\ND.F.LE.(N  +  .5))  DErD=l. 

IF(F.GT.(N  +  .5)AND.FXE.(N+ 1.))  DETD^  1.+  (M-1.)/NZ*(IZ-(2.*N+ 1.)* 

P  NZ) 

EiETURN 

C  0<=Z0<=L  Z=Z1<Z0 
21  F=-IZ/2./NZ 

N  =  INT(F+ .49999) 

T=0. 

IH((N*.aLLOG(M)).LT.(10.‘ALOG(lO.)))  T=  iy(M**(2*N)) 
lF{N.EQ.O)  T=l. 

IF(F.GT.(N-.5)J\ND.F.LE.N)  DErD=T 

n'(F.GT.NAND.F.LE.(N  +  .5))  DETD=T‘(1.+  (M-1.)/M*(IZ+2.‘N*NZ)/NZ) 

RETURN 

END 


5.  CONCLUSIONS 


The  initial  report  by  Anan'ev  was  the  'geometric'  view.  This  view  holds 
that  the  unstable  resonator  has  a  ray  that  has  the  required  integer  number 
of  waves  in  round  trip  at  locations  outside  the  optical  axis.  The  path  is 
longer,  so  the  frequency  shift  is  down  (to  longer  wavelengths).  The  rays 
are  not  normal  to  the  surface  of  the  mirrors,  but  some  feedback  is 
expected  due  to  diffraction  of  the  small  spot.  Al  Paxtone  at  PL/LIDA  had 
espoused  such  a  view.  The  preferred  ray  is  one  that  splits  the  angle 
equally  at  each  side  of  the  resonator  and  has  the  required  frequency  to 
give  integer  waves  of  optical  path  in  a  round  trip.  The  spontaneous 
emission  meeting  this  criterion  will  be  retained  in  the  resonator  and 
amplified  to  much  higher  intensity  than  that  which  does  not  meet  this 
criterion. 

The  problem  with  the  geometrical  view  is  that  it  does  not  require 
coherence  between  the  two  sides  of  the  resonator.  A  'modal'  approach, 
proposed  by  Greg  Dentes,  assumes  that  the  modes  are  coherent  with 
respect  to  themselves  on  the  opposite  sides  of  the  resonator.  This  view 
holds  that  the  whole  mode  can  be  fed  by  spontaneous  emission  all  over  the 
resonator  into  its  adjoint  mode  (converging  wave).  The  propagation  of  the 
beam,  which  keeps  clear  fringe  patterns  as  it  goes  from  near-field  to 
far-field,  suggests  that  the  modal  view  is  superior.  Both  views  are 
compatible  since  the  modal  view  predicts  hot  spots  exactly  where  the 
geometric  view  finds  its  resonant  rays.  The  difference  is  that  the 
geometric  view  holds  the  two  spots  uncorrelated,  where  the  modal  view 
holds  them  to  consist  of  two  modes,  one  in  phase  and  one  out  of  phase. 
Spontaneous  emission  and  gain  at  the  core  of  the  resonator  should  make 
the  in-phase  mode  dominant.  Professor  Anan’ev’s  March  1995  report 
presents  a  multiple  pass  geometric  viewpoint  that  bridges  the  modal  and 
geometric  views. 

An  interesting  variation  of  the  modal  view  comes  from  comparing  the 
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situation  to  ordinary  dielectric  waveguides.  This  is  most  appropriate  for 
the  situation  where  the  divergence  in  the  resonator  is  produced 
continuously  by  a  transverse  index  of  refraction  variation  or  an  imbedded 
lens  train,  rather  than  by  curved  reflectors  at  the  ends  of  the  resonator. 
The  difference  is  that  the  modes  are  evanescent  in  the  core  of  the 
resonator  and  radiate  in  the  outer  regions.  This  is  just  the  opposite  of 
the  situation  in  ordinary  dielectric  waveguide  modes. 

The  significance  of  these  points  occurs  when  the  laser  is  switched  from 
below  threshold  to  above  threshold  operation.  The  existence  of  a  mode 
inside  the  cavity  might  effect  the  speed  of  the  turn  on,  and  a  structured 
pattern  might  seed  filaments  in  a  diode  laser  medium.  The  sum  of  all  the 
frequencies  is,  however,  fairly  smooth,  so  the  gain  saturation  due  to  the 
sum  of  the  modes  does  not  itself  produce  filamentation.  The  laser  modes 
that  do  turn  on  as  the  gain  is  increased  lie  at  the  central  axis  of  the 
herringbone  pattern.  These  modes  grow  to  fill  the  lateral  extent  of  the 
resonator,  replacing  the  light  in  the  below  threshold  spontaneous  emission 
driven  modes.  Spontaneous  emission  that  does  not  fall  into  the  preferred 
frequencies  is  quickly  rejected  by  the  resonator.  Thus  the  existence  of 
these  patterns  does  not  significantly  affect  the  performance  of  a  laser 
above  threshold.  Their  utility  may  be  as  a  diagnostic  tool  to  evaluate  the 
quality  of  the  processing  used  to  manufacture  these  lasers. 
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Introduction  &:  Overview 

-  Geometric  motivation  for  converging  and  diverging  modes 
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Unstable  Resonators 


Basic  Design  Idea  —  Geometric  Modes 
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Spectra!  Pattern 

InGaAs/AlGaAs  CW  URSL  200  X  500  pm 


@  0.75  X  Ith 


9670  A  @  0.9  X  Ith 
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500  um  Long  URSL,  Mag=2.5.  Mode  Spacing=2.4  Angstrom 
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Spatio-speclral  pattern  of  a  cold  cavity  unstable  resonator  (200  x 
500  |xm) 


HALF  SYMMETRIC  UR  :  R=2200um  L=500um 
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Neumann  series  solution  method: 
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Fixed  ga 


ASE  Patterns 
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Grey  Scale  Combinations 
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case:  S  =  Uq  (diverging  mode  input)  8c  Power  gain  on  mode 


S  =  optimum  (converging  wave)  input 
Power  gain  on  mode  U,,  =  G 
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K  (Ai^cc  )  K  =  Petermann 
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